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Abstract 

The interaction between spherical magnetic nanoparticles is investigated from micromagnetic 
simulations and ananlysed in terms of the leading dipolar interaction energy between magnetic 
dipoles. We focus mainly on the case where the particles present a vortex structure. In a first step 
the local magnetic structure in the isolated particle is revisited. For particles bearing a uniaxial 
magnetocrystaline anisotropy, it is shown that the vortex core orientation relative to the easy axis 
depends on both the particle size and the anisotropy constant. When the particles magnetization 
present a vortex structure, it is shown that the polarization of the particles by the dipolar field of 
the other one must be taken into account in the interaction. An analytic form is deduced for the 
interaction which involves the vortex core magnetization and the magnetic susceptibility which are 
obtained from the magnetic properties of the isolated particle. 
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1 Introduction 



With the increasing progress in the synthethis of magnetic objects of nanometric scale such as spherical 
nanoparticles, nanodots, nanorings or layered films the diversity of systems made of such nano-objects 
as building blocks either as 2D or 3D assemblies in non magnetic environment [1-4] or in colloidal 
suspensions as ferrofluids [5,6] is continously growing. The magnetic behavior of magnetic nanometric 
particles either isolated or in nanostructured bulk materials is now quite well undertood both from 
experiments or numerical calculations [1,3,7] but a precise knowledge of interparticles interactions and 
of their influence on the macrocopic properties is still needed. Indeed, the interparticle coupling has 
been investigated in a variety of systems, such as nanograins [8-12] nanorings [13-15] or cylindrical 
nanodots [16-19] with a predominant attention paid on short range effects, such as exchange coupling, 
or on the influence of the coupling between single domain particles on the global magnetic properties. 
Conversely the long ranged interaction still deserves attention especially in cases where the magnetic 
structure of the isolated particle is complex (vortex [20-24] or onion [25,26] states for examples). A lot of 
work remains to be done on this point especially for spherical particles; in particular it seems important 
to develop models including the long ranged and anisotropic dipolar interaction. In the simple case of 
single domain particles the leading term in the interaction is the long range dipolar interaction which 
may lead to complex structures according to the shape of the particles on the one hand and the density 
and the dimensionality of the whole sytstem on the other hand [27-30]. In the case of particles with 
a non trivial internal magnetic structure, the interaction between particles is to be determined first. 
Indeed, it is generally admitted that when particles present a vortex structure, the resulting strong 
reduction of the magnetic moment at zero external field makes the dipolar interaction negligible. One 
aim of this work is to examine this point more precisely. In this work we focus on the interaction 
between spherical particles made of soft magnetic material (permalloy as an example) when they reach 
the vortex regime. We consider the simple situation of only two approaching spheres in a dumbell 
configuration. We especially compare the caculated interaction to the expected dipolar term. It is 
shown that two parameters characterizing the isolated particle play a central role: the magnetization 
of the vortex core and the suceptibility from which the polarization energy of one sphere in the dipolar 
field due the second one is calculated. 



2 Magnetization structure and hysteresis 

Since our purpose is to model soft magnetic particles in a general way rather than to focus on particles 
of a given material, the magnetic characteristics are somewhat arbitrary and correspond roughly to 
permalloy: the value of the exchange constant is A x = 1.10 J/m, the saturation magnetization 
J s = IT and the anisotropy is of uniaxial symmetry with a constant K\ ranging from K\ = to 
K\ = 7.10 4 J/m 3 . The particle radius is varied from R = 10 nm to 45 nm. In the following the particle 
volume will be denoted by v s . We determine the magnetization structure in the framework of the 
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micromagnetism formalism from the minimization of the total energy which is given by 



Etot — E x + E a + Edm + Ez 

A x Y,i(X7mi(r)) 2 + Ki(l - (m(r).a) 2 ) - ^p M s m(r).H dem - p M s H ex .m(r) dr (1) 

where E x , E a , E dm and Ez are the exchange, anisotropy, demagnetizing and Zeeman terms respectively, 
a is the unit vector in the direction of the easy axis and Q is the total volume of magnetic material, wich 
can include more than one particle. m(r) is the reduced magnetization density, related to the total 
magnetization by Mt = Jq M s m{r)dr. The calculations are performed with the micromagnetic code 
MAGPAR [31] which is based upon a finite element method. The problem includes two length scales, 
namely the exchange length, l ex = (2(jlqA x /J% ) 1 / 2 and the Bloch wall length Ib = \J A x jK\. Here, 
given the parameters chosen we have l ex = 5.013nm and Is > Ylnm. The value of the dimensionless 
parameter, K = 2K\/(uqMs) as defined in [20] ranges in between and 0.175. The mesh used in 
the calculations is such that the largest tetaedron size is smaller than l ex which imposes typically a 
mesh with Nj e ~ 10 5 elements for one sphere of radius R ~ 40 nm. We first calculate the magnetic 
structure of one isolated sphere in terms of both the size and the anisotropy constant K\. In order 
to characterize the magnetic state in the particle, we consider the local magnetization profile, m(r), 
which in the vortex regime, is decomposed in its cylindrical components using the vortex axis, say v, 
as the cylindrical axis 

m(r) = m v v + m^ip + m p p (2) 

where p and <p are the radial and tangential unit vectors of the projection of r in the plane normal to 
v. In the following, hatted letters denote unit vectors. The axis v is defined and actually determined 
as the mean direction of the local magnetization in the central part of the vortex, as shown in figure 
([T]). In a first step we focus on the behavior of the magnetization M in terms of the external field, 
H ex , especially for the variation of the field from small values up to saturation field; however we do 
not focus on the nucleation field. First of all, as is well known, small particles up to a threshold value, 
Rsd, are uniformly magnetized in a single magnetic domain and the hysteresis curve is a square. With 
our set of parameters, we get Rsd = 18 nm, for K\ = and 22 nm for K\ = 3 10 4 J/m 3 , in agreement 
with the result of the micromagnetic calculations of [32] and with the estimation given in [20]. Then 
a vortex structure is obtained, characterized by a vanishing value of the radial component m p (r), and 
\ m f{p)\ varying from |m^,(p)| = inside the vortex core, p < r c , to \m v (p)\ = 1 in the vicinity of the 
particle surface p ~ p m ax = Rsin(9(z)). At zero external field, the vortex direction, v is arbitrary 
when K\ = 0, while for K\ ^ 0, the direction taken by v relative to the easy axis a is controlled by 
the anisotropy energy which tends to allign m(r) on a. The anisotropy energy depends on both the 
value of K\ and the volume fraction of the particle where rh{r) is oriented parallel or antiparallel to 
a (\m.d\ ~ 1). The ratio of the volume fraction corresponding to the vortex core v c , characterized by 
m oriented parallel to v, to the volume fraction where rh{r) is oriented normal to v is directly related 
to the volume of the particle, v s . Roughly speaking the ratio of the total volume where \m.a\ ~ 1 is 
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either v c /v s or (1/2) (v s — v c )/v s if v is parallel or normal to a respectively. According to this scheme 
the stability condition for the vortex direction v to be normal to a reads v c < (l/3)v s . We can refine 
this very crude determination of the thershold by introducing the magnetization profile and imposing 



in (J2J) m p = 0. Then = \J\ — m 2 , and we get 

*R r R(z) 



dz I m v {p) 2 2iipdp 
-R Jo 
R pR(z) j 



Et = -K x I dz I -(1 - m v (p) 2 )27T P dp (3) 
J-R JO 1 

for v parallel or normal to a respectively. Then we assume that the component m v (p) depends on p 
only through r* = p/r sc , r sc being the pertinent scaling length (either l ex for K\ = or a function 
of both l ex and Is otherwise) and we neglect its dependence with respect to z. Hence, exploiting 
m v (p > r c ) = 0, we set the upper bound in the integral over m v to oo and we write (E^ — E'J) in the 
form 



Ei - ElJ = -Kt (lv 8 - 3Rr 2 sc I^ with / = 

The stability condition for a vortex normal to a is now 

R 



m v (r*) 2 2irr*dr* (4) 




>\hz (5) 



which must be read as R > r sc ^/9I / (2ir) = R t h(Ki) when R is varied at constant K%, or conversely 
as r sc < RsJ2-k/ (91) when the role of the magnetic characteristics is investigated for a given particle 
size. Equation can be rewritten in a more convenient form for the practical calculatioJ/m 3 ns of 
I, namely : S* = v / 9P(r sc )/(2vr) < 1 where I* is given by Jl]) with r* replaced by (r/R) and the 
upper bound replaced by (r/R) max = 1. Of course I* is then dependent on the value of r sc which is 
emphasized by the notation I*(r sc ). In any case, such an estimation is not supposed to provide an 
accurate determination of the threshold value of either R or t sc for the orientation of v normal to d 
but to predict at a qualitative level the effect of either the particle size or the magnetic parameters 
on the direction taken by the vortex. We can scale the vortex radius, r c on the smallest of the two 
characteristic lengths, r sc = inf(l ex ,lB)', however, this scaling may be taken with care and instead we 
can consider, when Ib increases, a scaling radius in the form of a function r sc (l ex ,lB). Notice that when 
using S* instead of S one has not to explicit the dependence of r sc . Therefore, at both A x and J s kept 
constant, we deduce from equ. ((H) that v gets _L a when the sphere radius is increased at K\ constant 
or when K\ is increased at R constant. As we shall see in the following the stability condition for the 
orientation of v relative to a agrees with this qualitative conclusion. Notice that the orientation of the 
vortex relative to the axis of easy magnetization is also found to be size dependent in the case of the 
cubic anisotropy [32] : in this latter case, the vortex is parallel to the axis of easy magnetization in 
large spheres. 

The orientation of the vortex relative to a can be determined from the magnetization profile, m(f) 
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as well as from the magnetization curve in terms of the external field, M(H ex ) by chosing the direction 
of the external field, h ex either parallel or normal to a. Indeed, we expect the magnetization process 
to differ according to the direction of the external field relative to the vortex one. We keep in mind the 
well known behavior of the vortex in the flat cylindrical nanodots where the magnetization is found 
to result from the shift of the vortex core when the field is applied normal to the vortex direction. 
Here, in the case of v(H ex = 0) _!_ a we expect a similar behavior for small values of the external 
field when h ex = a, up to the rotation of the vortex core along the direction of the field for high 
values of H ex before the magnetization in the whole volume of the sphere becomes oriented along h ex . 
Morever, in this case, we expect to have no remanence in the direction of the field, since in the vicinity 
of H ex = there is no net magnetization normal to the vortex direction. On the other hand a non 
vanishing magnetization at H ex = 0, corresponding to the vortex polarization in the direction ilo = 
h ex will be obtained. Conversely, the magnetization curve corresponding to h ex _L a still for a sphere 
characterized by v{H ex =0) la will present the more usual shape of a loop located in between ±H C 
with a non zero remanence corresponding to the vortex polarization. 

Then we focus on the external field induced magnetization in the spherical particle. As is generally 
obtained in nanodots or spherical soft magnetic particles [33, 35, 36], the magnetization M in the 
direction of the external field is found to vary nearly linearly with respect to the field, at least in the 
vicinity of H ex = and of course away from switching points where the vortex reverses as a whole. 
Such a linear behavior is observed both when h ex — v or flex -L *v • (or equivalently h ex -L cl or hex — & 
when R > R t h{Ki)). This means that the susceptibility x defined as 

X (6) 



dH ex 

does not depend on the value of the field to a very good approximation. We emphasize that the 
suceptibility is well defined for particles in the vortex regime since no multidomain state occurs and 
therefore a demagnetized state at H ex = can be ruled out. Notice that the value of x depends on the 
direction of the field as will be discussed below, and we should distinguish \\\ from x± according to 
the direction of the field relative to v. When such a distinction is not necessary it will be omitted to 
lighten the notations and \ 1S to be understood as its value corresponding to the orientation chosen for 
the field. Only in the case of an external field direction h ex neither parallel nor normal to the vortex 
direction the consideration of the two values of % is necessary. The independence of % with respect to 
H ex can be exploited for obtaining the variation of the total energy with respect to the external field. 
We consider the variation of M starting from H ex = to a value of H ex such that no switching of 
the magnetization occurs up to H ex and we analyse the corresponding variation of the magnetization, 
AM as the polarization of the sphere induced by the field. We have AM(H ex ) = xH e x- On the other 
hand, we can deduce AM(H ex ) from the energy, E(H ex ) by writting an equilibrium equation 

dE tot (AM) = 

dAM 1 ' 

which determine the equilibrium value of AM. The total energy depends explicitly on H ex through 
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the Zeeman term, — noH ex (m(0)v .h ex + AM), where we have expressed the permanent magnetization 
in the absence of the field as M(H ex = 0) = m(0)v, m(0) being the magnitude of the vortex core 
magnetization in the absence of the field. Then from ([7]) we get 

d 

{Edm + E x + E a ) = (M)H ex (8) 



dAM 

Therefore we get the variation of the total energy in the form 



/■AM 

E = E(H ex = 0) + / (i H ex (AM')dAM' - ^H ex {m{Q)v.h ex + AM) (9) 
J o 

AM 2 

= E(H ex = 0) + lio— fi H ex (m(0)v.h ex + AM) (10) 

2x 



where we have used AM(H ex ) = x^ex- Notice that both ([8]) and (JQJ) are exact equations while ([TO]) 
holds only in the case of a linear dependence of AM(H ex ) with respect to H ex . The second term of 
the r.h.s. of J9)) or (fTO]) has a simple interpretation: it is the energy of polarization of the sphere and 
corresponds to the energy cost of the reorientation of the magnetization inside the sphere. Equ. fTOjl is 
to be compared to the expression of the energy density of an array of coupled dots presenting a vortex 
structure obtained in [34]; more precisely the polarization energy coincides with the second term of 
equ. (5) of Ref. [34] where the induced magnetization in the dot is related to the vortex shift, s. The 
polarization energy can be written equivalently as [io(AMH ex )/2 where AM is to be understood as 
the induced moment due to the external field. Finaly, x can be related to the variation of the energy 
minus the Zeeman term 

x= 1 d(E tot -E z ) 
VoH ex dH ex 

3 Interaction energy between magnetic spheres 

Now we focus on the determination of the interaction energy between two magnetic nanoparticles in 
terms of the interparticle distance, r\2. The interaction energy is defined in a usual way 

E int (l, 2) = E tot {l, 2) - E tot {r 12 oo) (12) 

where E to t denotes the total energy of the two particles system, and (1,2) is a short notation for the 
orientation and location variables of the particles when they are brought together. We expect to get 
a form dictated by the dipolar interaction between the magnetic moment of the approaching spheres 
which reads 

E dip = 47rr 3 — d 11 2(m 1 ,m 2 ,r 1 2) (13) 
d 112 (rhi, rh 2 , r u ) = rhx.m 2 - 3(mi.fi 2 )(m 2 .f i2 ) (14) 
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where rrtj are the magnitude of the magnetic moments, rhi and r\2 the unit vectors in the direction of 
both moments and of the vector joigning the two particles and d\\2 is the angular function characteristic 
of the dipolar intercation. In the case of single domain particles m, = M s v s where M s is the saturation 
magnetization of the particles, and the orientations m; result from the minimum of d\\2 given in JHJ). 
For particles without magnetocristalline anisotropy, this gives obviously : rh\ = m,2 = ^12 and dm = 
-2. On the other hand, if the magnetocristalline energy is non zero on both particles with easy axes di, 
the orientations rhi will result of the interplay between the anisotropy energy tending to align rhi on 
di and the energy (|13p tending to minimize the angular function (I14|) . Furthermore if K\ takes a non 
vanishing value only in one particle say i = 1 and if this value is large enough to impose rh\ = d\, rh,2 
must orient in the dipolar field due to particle 1 i.e. in such a way that du 2 = rh2-d\ — 3(m2.ri2)(ai.ri2) 
is minimum. 

Now we consider the case of particles large enough to present a vortex structure. In this case, the 
orientations of the effective moments of the particles are the vortex directions, and the values of 
the moments are no more equal to M s v s but correspond to the vortex cores magnetizations and can 
be obtained from the magnetization curves M(H ex ). Let us introduce the coefficients an = rrii/(M s v s ) 
(in the following, we shall only consider the case of identical particles, so we drop the index i). a 
depends of course on the location of the second particle, which will be denoted in short by a(r±2, dm) 
or a(l,2). The value taken by a is not trivial since on the one hand it must be determined from the 
characteristics of the isolated particle and from the polarization of the particle by the dipolar field of 
the second one. A simple approximation for the interaction energy can be built in the framework of the 
dipolar approximation by considering that each particle is in the dipolar field of the other one. Then, 
we have to take into account two contributions. The first one which corresponds to (fT3l) . is nothing 
but —miH(iip(ri2)(vi-hd,ip(2,l)) where H c iip(ri2)h c iip(2,l) is the dipolar field created at ri by particle 
at Y2 and the second one is twice the polarization energy of each sphere in the field of the second one. 
The second contribution has been introduced in (fTUj) for one particle in a constant external field. In 
the present case, the role of m(o) is played by M s v s a(oo) while the induced moment in the direction 
of the dipolar field is 



We first consider the case where the vortex Vi is free to orient in the direction of the dipolar field due 
to particle j ^ i. This is the most general case since it corresponds to both the absence of anisotropy 
or particle large enough for the vortex to be normal to the easy axis. In this case we have 



Now adding twice the second term of (fTUl) to the total dipolar energy we get for the interaction energy 



P = phdi P = xHdip(ri2)hdi P 



(15) 



phdip = (a( r i2,dm) - a{oo))M s v s h dip = Aa(r 12 ,du2)M s v s v 



(16) 



Eint(l, 2) 



fj, {M s v s ) 2 
47rr? 9 



a(oo)(a(oo) + Aa(r 12 , <2ii2))dii 2 (mi, m 2 , f\i) 



(17) 
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which coincides with the interaction energy between polar polarizable hard spheres [37]. This is the 
important result of this section. It relates the interaction energy to the magnetic charateristics of the 
isolated spheres, namely, a(oo) and x through Act. Notice that this form for the interaction energy 
should hold not only in the case of two particles but also more generaly for an assembly of particles. 
In the latter case, the solvation of the total dipolar field and thus the determination of Aa(ri2, dm) 
becomes a difficult task. In the simple case of two particles, introducing u = x/(47ri? 3 ) (= x*/3 where 
X* = x/ v s is the reduced susceptibility) we get 

Aa = — — ^ : — r (18) 

{(r 12 /R) 3 + ud U2 ) 

4 Results 

We analyse first the magnetic behavior of the isolated particle with a special attention paid on the 
characterization of the vortex structure at low external fields. The magnetization curve is displayed 
in figure ([2]) for K\ = and R = 45 nm and the corresponding magnetic structure, through the 
local magnetic moment components ([2]), is shown in figures [3] and [4] for the remanent state and in the 
vicinity of the coercive field, before and after the reversal of the vortex core. Since K\ = 0, the vortex 
direction, v coincides with the direction of the external field. These results put in evidence the vortex 
structure and in particular the vortex core is reversed as a whole at the coercive field, with a nearly 
frozen m(r) structure. Moreover, we find that the reversal of rh(f) results from a global rotation since 
the component m v changes sign. When K\ ^ as described at the qualitative level in section j2|) v is 
parallel to the easy direction a for small values of R, and becomes normal to a beyond a K\ dependent 
threshold value Rth(K\). When a _L v, the vortex core is free to rotate in the plane normal to a and 
therefore will orient parallel to the external field if h ex _L a. As an example, we show in figure j5]) the 
magnetization curve for the two directions of the external field h ex = a and h ex _L a, in the case R = 
45 nm. Moreover, in the former case, the magnetization parallel and normal to the external field, Mm 
and M± are displayed. The magnetization behavior in terms of the external field corresponds to the 
situation v _L a; indeed, the remanence vanishes when h ex \\ a, while Mj_ takes a nearly constant value 
when H ex is varied in the central part of the M±(H ex ) curve. Moreover this value coincides with the 
remanence obtained for h ex _L a or equivalently h ex \\ v and therefore corresponds to the vortex core 
magnetization. The independence of M± with respect to H ex in the central part of the M±(H ex ) curve 
shows that the variation of the magnetization M^(H ex ), i.e. in the direction normal to v : corresponds 
to a shift of the vortex core normal to the direction of the field. The vortex core magnetization is thus 
nearly constant and given by the value of M± in that part of the curve. This is in agreement with 
the magnetization process obtained in the flat nanodot vortex structures. Finally M± sharply vanishes 
when the vortex rotates in the direction of the field, where the magnetization curve Mm whith h ex = a 
presents the hysteretic wings, similar also to what is found in the flat nanodot case where however 
this last value of the field corresponds to the vortex anhihilation prior to the saturation of the dot. 
The behavior of M(H ex ) outlined above is coroborated by the evolution with the value of the field of 
the structure of m(r), shown on figure j6]), where we see that the linear variation of M»{H ex ) in the 
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central part of the curve can be associated to a shift of the vortex in a direction normal to rh. From 
the evolution of Mm for h ex = a with the particle size, displayed on figure J7)), one can determine the 
threshold value R t h{K\) beyond which v is normal to a. Here we find R t h — 28 nm for K\ = 3.10 4 J/m 3 
(strictly speaking, 26 nm < R t h < 30 nm). We have calculated numerically the integral I*(r sc ) defined 
after equation Jl} from which we find that the threshold condition © is satisfied (see table I| in 
good agreement with the onset of the vortex structure deduced from the magnetization. Indeed, from 
this calculation, we get S* = 1 for R = 26.5 nm when K\ = 3.10 4 J/m 3 and thus Rth = 26.5 nm 
in agreement with the value deduced from the behavior of the magnetization M(H ex ). Similarly by 
decreasing K\ at constant R = 45nm, we find that the range of external field where the vortex is 
normal to a is reduced and then vanishes for K\ = 2.10 3 J/m 3 . Therefore we confirm our prediction 
that v JL a for K\ > K\ t h(R) at constant R. Then the value of x is determined from the slope of the 
magnetization curve, M(H ex ) in terms of H ex . The results are listed in table [III We also check that 
equ. (fTTj) is satisfied (see table P). 



4.1 Interaction beween particles 

We first consider the case of monodomain particles; for this we chose R = 10 nm. As expected 
the interaction energy is exactly given by the dipolar term with rrii = M s v s . When K\ = for both 
particles, the energy minimization leads to dii2(l, 2) = -2 and we thus mainly test the l/r 3 2 dependence 
of the interaction. On the other hand, we have also considered the case where only one particle bears 
a non vanishing uniaxial anisotropy with a value of K\ large enough to impose the orientation of its 
moment, m parralel to the easy axis a. Then the moment of the second particle orient itself in the 
field of the fixed particle in order to minimise the angular function du2- This provide an additionnal 
test of the behavior of the interaction through its angular dependence. The results are displayed in 
table IIIIl Now we consider the vortex regime with particles of radius R = 35nm or R = 45 nm. We 
start from particles without anisotropy, K\ = 0. In this case only one value for the susceptibility, X||, 
is to be considered, since the vortex allign spontaneously in the direction of the dipolar field. The two 
parameters involved in the expression of the interaction, x an d ce(oo), are determined first from the 
magnetization curve of the isolated particle. As a first test, we look at the angular dependence of the 
interaction energy. To this aim we start from the two spheres at a large distance and we minimize the 
total energy corresponding to non interacting spheres. Then we decrease the distance ri2 down to a 
not too small value of the ratio ryi/R and we perform a rotation of one sphere, say 2, arround the 
other one which is kept fixed. In this first calculation, we just calculate the components of the energy 
without minimization; we thus obtain the energy at a fixed value of local magnetic structure in the 
spheres, disregarding the polarization energy. The result is displayed in figure (jSJ) in the case r\i/R = 
4 and different values of the angular function d\\2 calculated by using m, = t),. We clearly obtain a 
linear dependence of Ei n t(l, 2) in terms of ^112(1, 2), and moreover the proportionality factor is exactly 
the result of the dipolar interaction, as deduced from (TT3]) . We thus conclude that when the structure 
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inside the spheres is frozen, the resulting interaction energy is indeed given by the dipolar interaction 
between the vortex cores. Then we consider the interaction energy after relaxation of the structure in 
the spheres, namely from the result of the total energy minimization in terms of the distance between 
particles. As expected and in agreement with eq. (PTT1) the value we get for dn2(l,2) is very close 
to d\i2 = -2 especially for short distances. The result for the interaction is shown in figure ((9|). We 
also compare the result corresponding to the dipolar interaction including the polarization energy or 
without this las term. This later approximation amounts to model the interaction by that between 
the dipoles corresponding to the isolated particles vortex cores. The approximation introduced in (TT71) 
is in very good agreement with the calculated result, for distances down to rn/R ~ 2.75, and the 
agreement for r^/R = 2.5 is still fairly good. Moreover we see that the inclusion of the polarization 
energy is quite important; indeed, the dipolar interaction calculated with the moments resulting from 
the isolated particles vortex cores reproduces the interaction only for distances larger than 3.35-R. 

In the case of particles with non zero uniaxial anisotropy, we focus on a situation where the vortex 
direction, v, is normal to the easy axis at zero external field. As an example we choose K\ = 3.10 4 J/m 3 
and either R = 45 nm or R = 35 nm. One can impose the plane in which the vortex is free to rotate 
via the direction chosen for the easy axis. Here we consider two situations where the two particles 
have the same easy axis, say a = z and the unit vector joigning the particles fi2 is either normal or 
parallel to a. Thus the equilibrium configuration of the particles corresponds to v\ = V2 = fi2 and 
dii2(l,2) = -2 in the former case and v\ = - Vz -L fi2 and dn2(l,2) = -1 in the latter case. The 
results are summarized in figure (fTUl) where we plot the interaction energy normalized by the value at 
the shortest distance considered, r\i = 2.25-R. The interaction energy is still very close to the dipolar 
plus polarization energy, eq (fl~7l) when d\vi = — 1, while in the case where the vortices are in line, the 
agreement for short distances is more qualitative. This is mainly due to an underestimation of the 
induced polarization by the dipolar field. We are lead to this conclusion by fitting the values of the 
parameters a and u in order to reproduce the calculated interaction energy by equ. (|17p . Doing this 
we can reproduce the calculated interaction energy only by using a non negligible enhancement of u 
while the fitted value of a remains very close to that calculated on the isolated sphere. The fitted 
results are also displayed in figure (fTUl) . To get a similar agreement with what is obtained in the case 
#1 = with d 112 = -2, the fitted value of u and a are 1.25u ca/c , 1.03a ca/c and 1.45u ca/c , 1.05a calc for 
R = 35 nm and 45 nm respectively. 

The results of this work are twofold. First we have precised the local magnetic structure in the 
sphere, and shown that beyond the well documented single domain to vortex transition in the case 
of a uniaxial anisotropy the vortex direction is normal to the easy axis once the particle radius is 
larger than a threshold value, R t h(Ki) for which a simple estimation is given. Then the interaction 
between particles is shown to present a dipolar character depending on two parameters characterising 
the isolated particle, namely the vortex core magnetization and the suceptibility. The vortex core mag- 
netization is strongly reduced when compared to the saturation magnetization M s which is quantified 
by the parameter a ~ 0.2 and this makes the interaction rather small but nevertheless non negligible. 
The order of magnitude of the interaction energy at distance r\<i = 2.25i? is slighltly smaller than the 
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barrier necessary to reverse the vortex core. However due to both its long range and its anisotropy the 
dipolar interaction is likely to lead to measurable effects in experimental assemblies of such particles. 
On a qualitative point of view, we do think that some of the finding of [24] are in agreement with the 
manifestation of dipolar effects, namely the tendency to form chains and to allign the vortex cores. In 
ref ( [24]) a micromagnetic simulation was already performed and was in agreement with the experi- 
ments; however, here we go a step forward by clearly pointing the dipolar character of the interaction 
between spherical nanoparticles. This allows us to predict that in a general way the behavior of dipo- 
lar and polarizable hard spheres will be transferable to assemblies of such particles even in the vortex 
regime. In this field, a very rich panel of structures is expected for both 2D systems [27] (and reference 
therein), [28,29] and 3D systems [30]. 
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Table I: Value of S* = v / 9/*(r sc /(2vr) involved in the stability condition ©. K x = 3. 10 4 J/?n 3 . J* 
is defined by equ. J4]) with r sc = R. 



R (nm) 


26 


30 


37 


45 


S* 


1.042 


0.685 


0.596 


0.533 



Table II: Reduced magnetic susceptibility calculated from (a) : equation ([6]); (b) equation (TTTj) - 





R 




x * (6) 









45 


3.229 


3.21 









40 


3.288 


3.28 









37 


3.336 


3.33 









35 


3.352 


3.37 






3.10 4 


45 


2.887 


2.94 


4.589 


4.58 


3.10 4 


35 


2.990 




6.510 


6.676 



Table III: Angular dependence of the interaction between modomain particles. R = 10 nm; -f^i(l) = 
7.10 5 J/m 3 i^i(2) = 0. @i and ©2 denote the angles (ai, i), (mi, z) and (m2, -2) respectively. 

r i2 = 4i2. dx™™^ is the minimum value of the angular function d\\2 corresponding to Gi fixed and d^^ 1 
is the result of the numeriacal calculation. Ei nt is the interaction energy per unit volume. According 
to the dipolar interaction the theoretical value for Ei nt /d\\2 is 2072 J/m 3 . 



0(ai) 


0! 


e 2 


,(caZc) 
u 112 


,(min) 
"112 


E int /dii2 (J/m 3 ) 


vr/8 


vr/7.948 


vr/4.442 


-1.2019 


-1.202 


2012 


vr/4 


vr/3.987 


vr/2.888 


-1.5828 


-1.584 


2007 


ir/2 


vr/2 


vr/2 


-2.0 


-2.0 


2004 
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Figure captions 

Figure [1] Local magnetization structure in the vortex regime. (R = 45 nm; K\ = 0; remanent state.). Top 
: projection of the local magnetic moment in the equatorial plane of the sphere, normal to vortex 
axis, v = M/\\M\\. The length of the arrows is proportional to the norm of the projection of m(r), 
rhpif). The central part of the vortex is clearly identified as the region where rn p (r) = 0; the 
direction of magnetization in this region coincides with the vortex direction, v. Bottom : Local 
magnetization in the direction normal the equatorial plane of the vortex shown on top, along a 
diagonal of this last one. The vortex direction is shown as the large bold arrow. 

Figure [2] Magnetization curve in the direction of the field. R = 45 nm; K\ = 0. 

Figure [3] Cylindrical components of the local magnetization profile at remanence. across the equatorial 
plane of the sphere (z = 0). m v , solid line; m^ dashed line; m p dotted line. R = 45 nm, K\ = 0. 

Figure H] Components m v (triangles), m v (squares) and m p (circles) of the local magnetization profile in 
the vicinity of the coercive field before (solid symbols) and after (open symbols) reversal of the 
vortex core. R = 45 nm, Ki = 0. d = p sign(y) where p is radius in the equatorial plane (z =0). 

Figure [5] Magnetization curve parallel and normal to the field . h ex = a (solid line) or v (dashed line). 
R = 45nm, K x = 3.10 4 J/m. 

Figure [6] Components m v (solid line), m y (dashed line) and m p (dotted line) of the local magnetization 
profile relative to the vortex core across the (x = 0) plane. R = 45 nm, K\ = 3 10 4 J/m 3 , H ex = 
68 kA/m, a = h ex = z, and v = x . The vortex is shifted along the y axis in the y > direction, 
by an amount y c = 15.75 nm, leading to a non symmetric range of variation for d. The location 
of the vortex core is indicated by the arrow. 

Figured Magnetization curve parallel and normal to the field for Ki = 310 4 J/m 3 and different sizes. 

Magnetization in the direction of the field and : R = 26 nm (solid); 30 nm (short dash); 37 nm 
(long dash). Magnetization normal to the field : R = 30 nm (dot short dash); 37 nm (dot long 
dash). For R = 26 nm, the magnetization normal to the field vanishes and the magnetization 
reversal occurs at a positive field since R = 26 nm enters in the range of particle sizes where the 
vortex direction is parallel to the easy axis, chosen as the direction for the field. 

Figure [8] Variation of the interaction energy for 2 spheres at r\2 = 4i? with the angular function d\\2 
characterizing the relative orientations, normalized by its maximum value, E{d\\2 = 2) — E(0). 
R = 35 nm; K\ = 0. 

Figure [9] Interaction energy per unit volume between two approaching spheres. R = 35nm; K\ = 0; 

d\\2 = —2. Open triangles: full calculation, from (fT2|) (the thin line is a guide to the eye); solid 
line: equ. (|17p : dashed line: simple dipolar approximation, u = 0. 

Figure [10] Same as figure j9]) for the interaction normalized by the value at r\2 = 2.25R. dotted lines: result 
of equ. ffTTl) with the values of a and u fitted in order to improve the agreement with simulated 
results. Open triangles : K\ — 0, d\\2 = -2, R = 35 nm; open squares : K\ = 3 10 4 J/m 3 , d\\2 
= -2, R = 35 nm; solid triangles : K\ = 3 10 4 J/m 3 , d\\2 = -1, R = 45 nm; solid squares : K\ = 
3 10 4 J/m 3 , d\\2 = -2, R = 45 nm. The different curves are shifted along the r\2 axis for clarity. 
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Figure 8: 
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Figure 9: 
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